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Abstract 


In this paper, dynamic behavior and performance of a fuel cell power plant (FCPP) which operates in parallel with a battery bank is tested 
under classified load conditions, such as mostly resistive, mostly inductive, resistive-inductive and non-linear loads. Thereafter, voltage stability 
analysis is performed using the dynamic response of the FCPP for stand-alone residential applications. Simulation results are obtained using the 
MATLAB® and Simulink® software packages, based on the mathematical and dynamic electrical models of the system. Using the experimental 
results, a validated model has been realized and voltage stability analysis is performed through this model. 
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1. Introduction 


Fuel cell power plants (FCPPs) are electrochemical devices 
that convert the chemical energy of a reaction directly into the 
electrical energy. Among the various next generation power 
plants, the proton exchange membrane fuel cell (PEMFC) sys- 
tems has been found to be one of the most promising energy 
sources due to high efficiency and environment friendly opera- 
tion [1]. It is important to understand and evaluate the dynamic 
response of PEMFC systems to current, voltage and power 
changes for controlling, designing, and optimum operation of 
the system [2]. As reported in Ref. [3], a FCPP alone may not be 
sufficient to satisfy all the load demands, especially during peak 
demand periods or transient events. Over the past few years, 
various studies performed to understand the dynamic behavior 
of PEMFCs. For example, Reference [4] discusses the behavior 
and performance of a PEMFC stack under fast load commuta- 
tions, and Ref [5] describes a small signal model of the PEMFC. 
Various studies have been reported in the literature with regard 
to the dynamic behavior of the PEM fuel cell power plants under 
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various loading conditions. Among them, Benziger et al. inves- 
tigated the dynamic response of PEM fuel cells to changes in 
load from a chemical point of view [6]. In Ref. [7], the cur- 
rent distribution and temperature distribution as well as their 
dynamic changes in fuel cell stack have been evaluated in situ. 
The information on the internal behavior of three small polymer 
electrolyte membrane fuel cells under static and dynamic load 
conditions was extracted in Ref. [8]. Yan et al. investigated the 
steady state and dynamic performance of PEM fuel cells under 
various operating conditions and load changes in Ref. [9] using 
the steady-state polarization curves, transient J-V response and 
electrochemical impedance spectroscopy techniques. However, 
there is a lack of experimental investigations have been per- 
formed to study the dynamic behavior of stand alone FCPPs 
under sudden load switching and related load conditions such 
as active-reactive power loading and loads with high harmonic 
distortion. 

It is important to understand FCPP’s dynamic behavior in 
order to estimate the time constants and to understand how 
a FCPP system with a reformer accommodates transient load 
conditions. Besides, it is also important to get an insight how 
the reactive power is supplied in FCPPs and how the inductive 
loads influence the system. Furthermore, we are also interested 
to know the dynamic behavior of FCPPs under load conditions 
consisting of harmonic components. 
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As FCPPs are expected to be widely used in the near future, 
voltage stability analysis of these next generation power systems 
must be considered during planning, operation, and control sce- 
narios. This will ensure that FCPPs will supply more stable, 
reliable, and high quality energy. Traditional voltage stability 
analysis can be modified and adapted to the FCPP based power 
systems. 

Voltage stability is associated with the capability of a power 
system to maintain acceptable voltages under normal and 
adverse conditions. On the other hand, voltage stability is also 
defined as the ability of the system to provide adequate reactive 
power support under both steady-state and transient conditions 
[10,11]. In voltage stability analysis, load characteristics play an 
important role [12], and both active and reactive power demand 
affects system voltage stability [13]. Also, active and reactive 
power losses increase when an increase in reactive power trans- 
fer occurs [13]. 

To find stability limits via voltage stability analysis, Newton— 
Raphson approach for power flow analysis may be used. This 
stability limit, also known as critical point, is often designed as 
the point where the power flow Jacobian is singular [14] i.e., 


det[J] = 0. () 


A 5kW FCPP equipped with a battery bank and a power con- 
ditioning unit (PCU) is currently operated in authors’ laboratory. 
In this paper, this FCPP is tested under various load conditions in 
order to obtain the dynamic response of the system. The classi- 
fied loads are switched for predefined durations when the FCPP 
is in stand-alone operation mode and the dynamic behavior of 
the system is observed. Using the experimental results, a val- 
idated model has been realized and voltage stability analysis 
is performed through this model. While experimental results 
are obtained using Fluke® power analyzers and PlugPower® 
FCPP system interface, simulation results are obtained using 
the MATLAB and Simulink software packages. 


2. System description and methodology 
2.1. Voltage stability analysis 


When the Jacobian of Newton—Raphson power flow equa- 
tions becomes singular, the voltage stability limit i.e., critical 
point of the system can be determined [15]. The parameters 
used in the voltage stability analysis are as follows: 


Vac ac output voltage of the inverter [V] 
VL load terminal voltage [V] 

PL load active power [W] 

QL load reactive power [VAr] 

Scritical Critical angle [rad] 


Vi -critical Critical load voltage [V] 
PL -critical critical load power [W] 


X reactance of the line connecting the combined system 
[2] 
ô phase angle between Vac and V [rad] 


7) load phase angle [rad] 


The active and reactive power transfer equations can be used 
to determine the critical parameters for voltage stability analysis 


Vac Vi 
po 2 L si 
X 


n ô, (2) 
v — Vac VL cos ô 


QL = 7 


(3) 


The critical values are determined by considering the singu- 
larity of the Jacobian matrix using the Newton—Raphson power 
flow analysis for the two-bus system, expressed as 


Vac VI 
fi(Vac, VL, ô) = PL — e- sin ô, (4) 


Ve — Vac V cos ô 


f2(Vacs VL, 8) = QL E 


(5) 


Eqs. (4) and (5) can be rearranged to form the Jacobian matrix 
as 


afi əfi 
AP | | a av |] Aê D 
AQL| |32 af2 | | AVI” 

ə OV, 


The singularity of the Jacobian matrix may be obtained by 
setting the determinant of Eq. (6) equal to zero i.e., 


of a ah ð 
fi f2 f2 afi = (7) 
06 OV, 06 ƏVL 


and, if we put Eqs. (4) and (5) into Eq. (7) then 


Vac = 2V1 cos ô. (8) 


Thus, we get the relationship between inverter output and 
load bus voltages at the critical point. To determine the critical 
angle, the dependence of Vac and ô is used instead of VL. The 
relationship between the real and reactive powers, PL and QL 
may be expressed as 


QL = Pi tan Q (9) 
Using Eqs. (2), (3) and (9), the critical angle may be expressed 
as 
—1 


tan (28) = ene 


(10) 


Finally, the critical angle and power values can be obtained 
from Eqs. (2), (8) (9) and (10) as 


1 —1 
Ôcritical = 7 arctan (; ) (11) 


an Q 
and 


2 
Vie tan Scritical 


2X f aa 


PL-critical = 
respectively. The phase angle, whose critical value is obtained 
earlier, depends on the amount of power that the FCPP produces 
and on some other parameters explained in Section 2.2. 
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2.2. Mathematical modeling of PEMFC system 


The mathematical model parameters of the FC system used 
in this model are as follows: 
B,C constants to simulate the activation over voltage in 
PEMFC system and [V] 
CV conversion factor [kmol of hydrogen per kmol of 


methane] 
E Nernst instantaneous voltage [V] 
Eo standard no load voltage [V] 
F Faraday’s constant [C kmol~!] 
Tec FC system feedback current [A] 
kı PI gain 
Kan anode valve constant [,/kmolKg(atm s)~!] 
Kp, hydrogen valve molar constant [kmol (atm s)] 
Kmo water valve molar constant [kmol (atm s)!] 
Ko, oxygen valve molar constant [kmol (atm s)” 1] 


K, modeling constant [kmol (s A)7}] 
Mp, molar mass of hydrogen [kg kmol" !] 


No number of series fuel cells in the stack 
PH hydrogen partial pressure [atm] 

Pmo Water partial pressure [atm] 

POs oxygen partial pressure [atm] 

qHz molar flow of hydrogen [kmol s~!] 
qo input molar flow of oxygen [kmol s~!] 


methane methane flow rate [kmol s~!] 

di, hydrogen input flow [kmol s~!] 

qi hydrogen output flow [kmol s~!] 
IH, hydrogen flow that reacts [kmol s~!] 


du amount of hydrogen flow required to meet the load 
change [kmol s7!] 

R universal gas constant [(1 atm) (kmol K)-!] 

Rint FC internal resistance [Q] 

T absolute temperature [K] 

U utilization rate 

Van volume of the anode [m°] 


Veell dc output voltage of FC system [V] 


T1, T2 reformer time constants [s] 

T3 time constant of the PI controller [s] 
TH hydrogen time constant [s] 

TO, oxygen time constant [s] 

THO water time constant [s] 

Nact activation over voltage [V] 

Nohmic Ohmic over voltage [V] 


The relationship between the molar flow of any gas (hydro- 
gen) through the valve and its partial pressure inside the channel 
can be expressed as [16] 


K 
ee ee ce (13) 


PH) V Mnp, T 


For hydrogen molar flow, three significant factors must be 
considered—hydrogen input flow, hydrogen output flow, and 
hydrogen flow during the reaction [16]. The relationship among 


these factors can be expressed as 


out 


RT in 
P= 7 (4H, — 4H, 
an 


ip). 14 
n da) (14) 
According to the basic electrochemical relationship between 
the hydrogen flow and the FC system current, the flow rate of 
reacted hydrogen is given by [16] 


r _ Nolrc 
M= oF 
Using Eqs. (13) and (15), and applying Laplace transform, 


the hydrogen partial pressure can be obtained in the s domain as 
[16] 


= 2K; Igc. (15) 


1/Ky, in 
= —— — 2K; lfc), 16 
aa ae mas He rl gC) (16) 
where 
Van 
TH, = : 17 
Ha = RT (17) 


Similarly, water partial pressure and oxygen partial pres- 
sure can be obtained. The polarization curve for the PEMFC 
is obtained from the sum of the Nernst’s voltage, the activation 
over voltage, and the ohmic over voltage. Assuming constant 
temperature and oxygen concentration, the FCPP output voltage 
may be expressed as [3,16,17] 


Vee = E + tact + Nohmic> (18) 
where 

Nact = —B ln (Clg) (19) 
and 

Nohmic = — R" Ipc. (20) 


Now, the Nernst’s instantaneous voltage may be expressed as 
[16] 


RT 
Eo + log 


N 2F 


PH V ra | , 21) 


The FCPP consumes hydrogen according to power demand 
while the reformer continuously generates hydrogen enriched 
fuel for consumption in the FC stack operation. The reformer 
model can be expressed as [16,17] 

qH, CV 


= A 22 
uns? +11 +Tr)s +1 (22) 


methanol 


During operational conditions, to control the hydrogen flow 
rate according to the output power of the FC system, a 
proportional-integral (PI) control system is used. To achieve this 
feedback control, FC current from the output is taken back to 
the input while converting the hydrogen into molar form [14,17], 
given by 


req _ Nolrc 
1m “SEU ` 


(23) 
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Table 1 

The parameters used for the mathematical modeling of the FC system 
Specification Value 

Activation voltage constant (B) 0.04777 [A7!] 
Activation voltage constant (C) 0.0136 [V] 


Conversion factor (CV) 2 

Faraday’s constant (F) 96484,600 [C kmol~!] 

FC system internal resistance (Rint) No x 0.00303 [Q] 

FC absolute temperature (7) 343 [K] 

Hydrogen time constant (ty, ) 3.37 [s] 

Hydrogen valve constant (Ky, ) 4.22 x 1075 [kmol (s atm)~!] 
Hydrogen—oxygen flow ratio (rH o) 1.168 

K, constant = No/4F 2.2802 x 10~7 [kmol (s A)~!] 


Line reactance (X) 0.0809 [Q] 

Methane reference signal (Qmethref) 0.000015 [kmol s~!] 
No load voltage (Eo) 0.8 [V] 

Number of cells (No) 88 

Number of stacks (Ns) 1 

Oxygen time constant (to, ) 6.74 [s] 


Oxygen valve constant, (Ko, ) 2.11 x 1075 [kmol (s atm)~!] 


PI gain constant (kı) 0.25 
Reformer time constant (T1, T2, T3) 2 [s] 
Universal gas constant (R) 8314.47 [J (kmol K)~!] 


Utilization factor (U) 0.8 
Water time constant (TH,0) 18.418 [s] 
Water valve constant (Ky,0) 7.716 x 1076 [kmol (s atm)~!] 


The amount of hydrogen available from the reformer can be 
used to control the methane flow rate by using a PI controller 
[16,17], expressed as 


ky Nolrc i 
qmethane = (x F =) ( = is) : (24) 


T3S 


The mathematical model of the FC system can be obtained 
using the above equations. The parameters for the mathematical 
modeling of the FC system are given in Table 1. 


2.3. Power conditioning unit (PCU) 


The FCPP produces dc electrical energy which must be con- 
verted to ac electrical energy for household appliances. The 
considered system for domestic loads requires a dc/ac inverter, 
which also acts as an energy management unit to match the dif- 
ferent dynamics of the FC system, battery bank and various types 
of loads. 

The parameters used in the mathematical modeling of the 
PCU in the mathematical model are as follows: 


Íbattery battery bank current [A] 

Tac de current [A] 

Ipc FC system current transferred to dc bus [A] 
I load current [A] 

m inverter modulation index 

Pac ac output power [W] 

PL load power [W] 

Vac ac output voltage of the inverter [V] 

Vac combined system dc terminal voltage [V] 
Vi reference voltage signal [V] 


VL load terminal voltage [V] 


X reactance of the line connecting the FC system and the 
load [Q] 
0 load phase angle [rad] 


The ac voltage and active power at the inverter output can be 
expressed as a function of the modulation index and phase angle 
as [16] 

Vac Vi 
Pyp = "sin (8). (25) 
X 

The corresponding load current and dc current can be 

expressed as 


PL 
IL =>: 
Vi cos (0) 
Tac = mI, cos (6 + ô), (27) 


(26) 


Assuming a lossless inverter, Pa, can be expressed as [16] 
Pac = Pac = Vac lac- (28) 


According to Eq. (23), the amount of required hydrogen flow 
can be calculated as 


reg NoUae + Ivattery) 
oe 2FU 
If the battery bank is discharging, the current is negative. If 
the battery bank is charging, the current is positive. According 


to Eqs. (2), (28) and (29) the variation of the phase angle 6 can 
be obtained as 


a. ((24u, FU/No) + Tpattery) X 
mV, : 


(29) 


ô= 


(30) 


2.4. Experimental setup 


The 5kW PEMFC system currently operated in our labora- 
tory feeds a 500-square-foot house as shown in Fig. 1. The FCPP 
can be operated in the stand-alone or in the grid parallel mode. 
In this study, experimental results are obtained when the FCPP 
was operated in the stand-alone mode. 

This FCPP’s PCU consists of an inverter, dc/dc converters 
and a battery bank. The PCU supplies power to internal auxiliary 
components of the FC system, manages charging and discharg- 
ing of the battery bank, and supplies the load requirements. There 
are four, 12-V maintenance-free lead-acid batteries integrated in 
series to the inverter. The battery bank supplies auxiliary power 
during cold-start of the FCPP, and supports peak demand peri- 
ods or transient events and compensates the tracking delays of 
the FC stack power. 

To observe the dynamic behavior of the FCPP and the effects 
of various loads on system voltage stability, the following types 
of loads were used: 


- Resistive loads, 

Inductive loads, 

- Resistive-inductive loads, 

Loads with high harmonic distortion. 
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Fig. 1. Test experimental setup: (a) FC system, (b) various loads and (c) measurement systems. 


The specific types of loads and the corresponding switching 
intervals are shown in detail in Table 2. Experimental results 
are obtained for the time interval between 04:02:20 P.M. and 
04:49:00 P.M. of the recorded data. This time interval corre- 
sponds to the 2800 s of duration. The chosen loads are switched 
according to the given order and all results are obtained accord- 
ing to the load types and time durations given in Table 2. 

The ac load bus quantities are measured using Fluke power 
analyzers while the quantities related to FCPP, PCU and bat- 
tery bank are measured using PlugPower FC system interface 
software. All data measurement and logging devices are syn- 


Table 2 
Various loads used in the experimental study 


Loads Time intervals (s) 


[70-259], [1134-1328] and 
[1338-1535] 

[266-332] and [339-466] 
[479-719], [950-1128] and 
[1338-1535] 

[735-1985] 


Washing machine 


Microwave 
4 x Computer (without monitors) 


4 x Small transformers (secondary 
open-circuited) 

1 x Small transformer (secondary 
open-circuited) 

2 x Small transformers (secondary 
open-circuited) 

Resistive load bank (switch 1) 

Resistive load bank (switch 2) 


[1841-1985] 
[1986-2105] 


[2113-2417] and [2608-2757] 
[2418-2607] 


chronized. At first, the FCPP was switched to stand-alone 
operation mode under a base load and the data logging proce- 
dure is started at t= 0. After 70 s, classified loads were switched 
on according to the sequence and durations given in Table 2. 


3. Results and discussions 


The block diagram of the integrated overall system is shown 
in Fig. 2. 

Figs. 3—8 show the results for the ac bus measurements involv- 
ing active power, reactive power, power factor, load voltage and 
total harmonic distortion (THD) of voltage and current, respec- 
tively. 

Figs. 3 and 4 show that there is a constant active and reactive 
power values of 130 W and 120 VAr because of the control and 
measurement devices connected to the ac side. The active and 
reactive power variations above these values are caused by the 
loads as given in Table 2. The load variations switched as shown 
in Table 2 was selected to build a variety of active and reactive 
loading combinations to test the PEMFC power plant under these 
conditions. Thus, we could observe the dynamic performance 
of the plant under different power factor values as shown in 
Figs. 3 and 4. According to the Figs. 3 and 4, resistive load 
bank draws mostly active power while transformers draw mostly 
reactive power. The washing machine and microwave draw both 
active and reactive power. 


M. Uzunoglu et al. / Journal of Power Sources 168 (2007) 240-250 245 


FC System 
Natural Gas 


Input f 7 5 2 


Auxiliary 
Power 


Power 
Conditioning 
Unit 


120V ac 


Charge 
& 


Discharge 


Battery 
Bank 


+ Collected, 


N S ., Sensor oe PlugPower System 
fe Pi, >, nterface Software 
me Bius da -o ee 
"es, > Data 

tee eeeeee > i 


Modem 


Single Phase 
Transformer 


ac 
Load Bus 


Washing 
Machine 


Microwave 


Computers 


Transformers 
(Inductive Loads) 


es Resistive 
«a T Load Bank 


Sensor 
Data 


Measurements 


Fig. 2. The overall system block diagram. 


From Fig. 5, we observe that the load bus voltage varies 
according to the load power which causes voltage drop on the 
equivalent system impedance between the inverter output and 
load buses. The reactive power oscillates between the inverter 
output and the load bus, thus reactive power demand can be 
met as a result of the voltage variations shown in Fig. 5. The ac 
load bus voltage oscillates between the 112.2 V minimum and 
116.8 V maximum. 


1800 r T r r 7 
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Q 
o 
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Load active power [W] 
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Time [s] 


Fig. 3. ac load bus active power variation under various load conditions. 


Nonlinear elements generally cause harmonics which result 
in current and voltage distortions in electrical systems. A mea- 
sure of the distortion is represented by the ratio of the equivalent 
harmonics to the percentage of the fundamental signal ampli- 
tude which is well known as THD. The THD of the voltage and 
current distortions are given in Figs. 6 and 7 respectively, when 
the loads were turned on in predefined time intervals. In partic- 
ular, the microwave and computers have the highest harmonic 
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œ 
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Fig. 4. ac load bus reactive power variation under various load conditions. 
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Fig. 5. ac load bus voltage variation under various load conditions. 
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Fig. 6. ac load bus total harmonic distortion for voltage variation under various 
load conditions. 
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Fig. 7. ac load bus total harmonic distortion for current variation under various 
load conditions. 
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Fig. 8. ac load bus power factor (PF) variation under various load conditions. 


distortions since they have power electronic circuits. The THD 
affects the power factor (PF) which considers the fundamental 
and harmonic components while traditional cos considers the 
fundamental component only. The variation of the PF is shown 
in Fig. 8. 

Figs. 9-15 show the results related to the FCPP, battery and 
PCU outputs, stack dc power, cell and stack voltages, battery 
power, inverter output voltage and current, average ac power 
and phase angle, respectively. 

To meet the required load power, FCPP produces the stack 
dc power shown in Fig. 9. This produced power includes the 
auxiliary power, charging power of the battery, and the power 
transferred to the dc bus. The FC stack produces dc power to 
track the active power demand of the load as depicted in Fig. 3. 
According to the Fig. 9, we can mention that the stack dc power 
follows the general trend of the active load demand. However, it 
does not match the demand at all because it has step-up and step- 
down time delays. The reasons of these mismatches are caused 
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Fig. 9. The variation of FC stack dc power according to load demand. 
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Fig. 10. The average cell voltage according to load demand. 
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Fig. 11. The FC system voltage according to load demand. 
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Fig. 12. Battery power including charges and discharges. 
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Fig. 15. The variation of ac voltage phase angle. 
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by the time constants of the fuel cell dynamics, particularly the 
time constants of the reformer dynamics. 

The response of the cell and stack voltages varies as shown in 
Figs. 10 and 11, respectively, according to the amount of power 
supplied by the FC system as shown in Fig. 9. The stack dc 
voltage and average cell voltage decrease when the FC output 
power increases due to the FC power—voltage relationship. The 
stack de voltage varies according to the active and reactive pow- 
ers drawn from the fuel cell. From Figs. 9-11, it is obvious that 
the PEMFC stack equipped with reformer has a slow dynamic 
response. 

The stack dc power given in Fig. 9 is insufficient to meet the 
load power requirements. It is obvious that the stack power only 
tracks the general trend of the power demand, and needs further 
conditioning to match the actual load demand. From Fig. 12, 
it is seen that the battery bank supplies the power to compen- 
sate all mismatches and tracking delays between the load and the 
stack powers by charging and discharging. Additionally, the bat- 
tery bank assists FC during peak demand periods and transient 
events. Without a storage device as an auxiliary power source, 
the FC power plant would have to supply all power demand thus 
increasing the size and cost of the FC power plant. The battery 
power becomes negative when it supplies power to the dc bus, 
and it becomes positive when it is recharged. The PCU uses 
stack and battery power to condition the output ac power. 

The inverter output voltage is successfully controlled by the 
PCU within the acceptable limits as illustrated in Fig. 13. How- 
ever, the most significant deviations of the ac output voltage 
(although it is not very high) is recorded when the loads includ- 
ing reactive components such as washing machine and small 
transformers are switched on. It may be mentioned that the mag- 
nitude of the voltage oscillations increases during the reactive 
load switching. However, the inverter keeps its output voltage 
within acceptable limits from minimum 119.5 V to maximum 
121.5 V while tracking the reference voltage. 

According to the inverter output active power (Fig. 3), reac- 
tive power (Fig. 4) and voltage (Fig. 13), the ac output current 
varies as given in Fig. 14. This current profile consists of active 
and reactive load currents. 

We also observe that this apparent current influences average 
cell and stack voltages given in Figs. 10 and 11, respectively. 
An increase in reactive power demand causes a proportional 
decrease in the FCPP output voltage. Although the fuel cell 
is a dc source, the resultant current of active and reactive cur- 
rent components has an effect on average dc cell voltage and 
stack de voltage which is an interesting result obtained in this 
experimental study. 

The variation of ô is depicted in Fig. 15 which can be attained 
using the parameters given in Eq. (30). The parameters of Eq. 
(30) are obtained using the experimental data and mathematical 
modeling. From Fig. 15, we can observe that the phase angle 6 
varies with respect to the active load demand from the fuel cell 
stack and it tracks the general trend of the active load variation. 

Figs. 16—19 show the results of critical angle and power at the 
load bus. The variations of critical values are produced exper- 
imentally through mathematical modeling to account for the 
parameters given in Eqs. (11) and (12). 
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Fig. 16. The variation of critical angle at the load bus. 
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Fig. 17. The variation of critical power at the load bus. 
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Fig. 18. The variation of critical angle at the load bus. 
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Fig. 19. The variation of critical power at the load bus. 


The variation of critical phase angle is given in Fig. 16. It is 
seen that the critical angle varies according to the power factor. 
All parameters which affect the power factor have an effect on 
the critical angle as well. Fig. 17 shows the variation of critical 
active power which is affected by the critical angle (ôcritica1), ac 
voltage and equivalent reactance of the line connecting the FCPP 
to the load. The critical angle and the critical power can also be 
called as the stability margins of the system. If we get closer to 
these critical values in any loading condition, the system will 
be very risky in terms of stable and safe operation. The voltage 
stability may be disturbed and voltage collapses may occur if 
these limits get enforced. 

Two of the above mentioned parameters are tested separately 
to ascertain how they affect the voltage stability. One of the 
parameters is the amount of reactive power drawn from the 
FCPP and the other is the reactance of the line connecting the 
FCPP to the load. At first, the reactive power is assumed to be 
two times greater then the original demand and then, the vari- 
ation of new critical angle is given in Fig. 18. Then the line 
reactance is assumed to be two times greater than the origi- 
nal value and the variation of new critical power is given in 
Fig. 19. From Figs. 18 and 19, it is evident that the increase 
in either the reactive power or the line reactance, results in a 
decrease in the critical power and critical angle. Thus, the sys- 
tem gets closer to the risky operating points in terms of voltage 
stability. 


4. Conclusions 


In this paper, the dynamic behavior of the FCPP is tested 
under various loads and their effects on voltage stability are ana- 
lyzed using a PEMFC based power system. Based on the analysis 
and experimental results, we can conclude the following: 


e Due to the slow dynamic response of the tested FC stack 
equipped with a reformer, FCPP needs an auxiliary energy 
storage system. The FCPP needs to be integrated with extra 


energy storage devices such as battery or ultra-capacitor bank 
to successfully meet all load demands. 

e Battery bank assists FCPPs during transient and peak demand 
periods, and tries to prevent the power tracking delays and 
mismatches. 

e To meet the reactive power demand results in voltage drop at 
the load bus. Both the active power and the reactive power 
demand influences the inverter output voltage, average cell 
voltage and stack voltage. 

e The power factor decreases due to the harmonic reactive pow- 
ers drawn by the nonlinear elements. Furthermore, it may 
be mentioned that the higher harmonic distortion negatively 
affects stability. 


From the voltage stability analysis, it is evident that the power 
and phase angle values do not reach the critical values in the cur- 
rent scenario. However, for different applications, the critical 
values must be considered in order to maintain the voltage sta- 
bility. For systems with long line distances and/or high reactive 
power requirements, the stability analysis should be performed 
prior to the design and operation of FCPPs. 
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